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Abstract 

Using data samples collected at ^/s = 3.686 GeV and 3.650 GeV by the BESII detector at 
the BEPC, the branching fraction of 'ip{2S) 3(7r'''7r“) is measured to be [4.83 ± 0.38 (stat) ± 

0.69 (syst)] X 10“^, and the relative branching fraction of J/'0 ^ 2(7r'''7r“) to that of J/if) lJ-~ 
is measured to be [5.86 ± 0.19 (stat) ± 0.39 (syst)]% via ip{2S) tt+tt” J/'0, d/b ^ 2(7r+7r“). 
The electromagnetic form factor of 3(7r'*“7r“) is determined to be 0.21 ± 0.02 and 0.20 ± 0.01 at 
a/s = 3.686 GeV and 3.650 GeV, respectively. 

PACS numbers: 13.25.Gv,13.40.Gp,14.40.Gx 
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I. INTRODUCTION 


Strong decays of to 3(7r’''7r“) are suppressed, since the reaction violates G-parity 

conservation. In colliding beam experiments, 3(7r’''7r“) may also be produced by 

> 7 * —3(7r’''7r“) (called the “continuum process” hereafter). It is expected that the 
continuum contribution is large and may contribute around 60% of the 3(7r’''7r“) events at 
the energy. This contribution must be removed in determining B{ip{2S) 3(7r+7r“)), 

as has been described for the 'ip{2S) tt+tt" decay mode [1]. 

In this analysis, data samples at the peak = 3.686 GeV) and off-resonance 

(y% = 3.650 GeV) are used. The continuum contribution at the peak is estimated 

using the off-resonance sample and subtracted to obtain a model independent measurement 
of the t/’(25') —> 3(7r’''7r“) branching fraction. We also use the samples to obtain the 3(7r’''7r“) 
electromagnetic form factor which allows us to calculate the branching fraction based on the 
theoretical assumption described in Ref. [1]. 

There is a big contribution from 'ip{2S) Jf-ip ^ 2 {tt^tt~) in our 'ip{2S) 

3(7r+7r“) sample. This process allows us to measure the branching fraction of J/ijj —> 
2(7r^7r“). The advantage of this method is that we need not subtract the continuum contri¬ 
bution for this process. 

The existing branching fraction measurement of 'ip{2S) 3(7r’''7r“) was done by the Mark- 

I experiment [2] based on (9±5) candidate events. The branching fraction oi J /'ip 2(7r’''7r“) 
was also measured by Mark-I [3] with (76 ± 9) events observed, and there is also a recent 
result for this decay reported by the BABAR experiment [4]. 


II. THE BES EXPERIMENT 

The data used for this analysis are taken with the updated Beijing Spectrometer (BESII) 
detector at the Beijing Electron-Positron Gollider (BEPG) storage ring. The 'ip{2S) data are 
taken at y/s = 3.686 GeV with a luminosity of £ 3.686 = (19.72±0.86) pb ^ [5] measured with 
large angle Bhabha events. The number of 'ip{2S) events is A^° 2 S) “ (14.0 ± 0.6) x 10® [6] 
as determined from inclusive hadrons. The continuum data are taken at y/s = 3.650 GeV, 
and the corresponding luminosity is £ 3,650 = (6.42 ± 0.24) pb ^ [5]. The ratio of the two 
luminosities is £ 3 . 686 /£ 3.650 = 3.07 ± 0.09. 
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The BESII detector is a conventional solenoidal magnet detector that is described in 
detail in Refs. [7, 8]. A 12-layer vertex chamber (VC) surrounding the beam pipe provides 
trigger and track information. A forty-layer main drift chamber (MDC), located radially 
outside the VC, provides trajectory and energy loss {dE/dx) information for charged tracks 
over 85% of the total solid angle. The momentum resolution is ap/p = 0.017a/1 -|- {p in 
GeV/c), and the dE/dx resolution for hadron tracks is ~ 8%. An array of 48 scintillation 
counters surrounding the MDC measures the time-of-flight (TOF) of charged tracks with 
a resolution of ~ 200 ps for hadrons. Radially outside the TOF system is a 12 radiation 
length, lead-gas barrel shower counter (BSC). This measures the energies of electrons and 
photons over ~ 80% of the total solid angle with an energy resolution of ue/E = 22%/ 

{E in GeV). Outside of the solenoidal coil, which provides a 0.4 Tesla magnetic held over 
the tracking volume, is an iron hux return that is instrumented with three double layers of 
counters that identify muons of momentum greater than 0.5 GeV/c. 

A GEANT3 based Monte Carlo (MC) program with detailed consideration of detector 
performance (such as dead electronic channels) is used to simulate the BESII detector. The 
consistency between data and Monte Carlo has been carefully checked in many high purity 
physics channels, and the agreement is quite reasonable [9]. 

In generating MC samples, initial state radiation is included, and 1/s or 1 /dependent 
form factors are assumed where required. MC samples of n^7i~ J/iIj, J/ijj X are 

generated with the correct 7r’''7r“ mass distribution [10], and J/tlj ^ 

is generated with the correct angle distribution. Other samples are generated 

according to phase space. 

III. MEASUREMENT OF tp(2S) 3(7r+7r-) 

A. Event Selection 

Six charged tracks with net charge zero are required. Each charged track, reconstructed 
using hits in the MDC, must have a good helix £t in order to ensure a correct error matrix 
in the kinematic fit. All six tracks are required to (1) originate from the beam intersection 
region, i.e. < 2 cm and |V;| < 20 cm, where 14 , Vy, and 14 are the x, y, and 2 : 

coordinates of the point of closest approach to the beam axis, and (2) have | cos 6^1 < 0.8, 
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FIG. 1: The invariant and recoil mass distributions of '4>{2S) 3(7r'''7r“) candidates at 

= 3.686 GeV. The events with two 7r''“7r“ pairs within (0.47,0.53) GeV/c^ are removed as 
KgKg background. While the events with a pair within (3.06,3.14) GeV/c^ are removed as 

■0(25) ^ 7r+7r“ J/0 background. 

where 6 is the polar angle of the track. 

A four constraint kinematic £t is performed with the six charged tracks assuming all of 
them to be pions. If the confidence level of the fit is greater than 1%, the event is categorized 
as 0(25') —3(7r’''7r“). 

Fig. 1 shows the invariant and recoil mass distributions of tt+tt" pairs. If the recoil 
mass of any 7r’''7r“ pair is between 3.06 and 3.14 GeV/c^, the event is considered a '0(25') — 
J/tlj, J/'^ 2(7r’''7r“) candidate and removed. If the masses of any two tt+tt" pairs is 

between 0.47 and 0.53 GeV/c^, the event is considered as '0(25') —^ n^'K~ 

background and removed. 

Applying these criteria to the data collected at ^/s = 3.686 GeV, 670 events survive, 
while for the data collected at ^/s = 3.650 GeV, 71 events remain. The efficiencies of these 
criteria are e^{ 2 S) = 6.8% for '0(25') —3(7r’''7r“) and Scont = 3.8% for 3(7r’''7r“). The 

lower Scant results from the initial state radiation correction (the maximum radiative photon 
energy is set to 0.7GeV) in the generator, which reduces the center-of-mass energy for many 
events generated. These events cannot survive the kinematic fit, which leads to the lower 
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FIG. 2; Confidence level distribution from the kinematic fitting, (a) is for ip{2S) 3(7r'*'7r“). The 
dots with error bars are data. The blank histogram is MC simulated signal plus the continuum 
contribution measured with ^/s = 3.650 GeV data (hatched histogram) and the MC simulated back¬ 
ground events (dark shaded histogram), after proper normalization, (b) is for e'’'e“ ^ 3(7r^7r“) 
measured with y/s = 3.650 GeV data. The dots with error bars are data, while the blank histogram 
is MC simulation. 


^cont • 

Remaining backgrounds in the t/’(2S') sample include: (1) residual KgKg or t/’(2S') — 
events; (2) events with kaons or electrons misidentified as pions, and (3) events 
with low energy neutral tracks like or 7 . Monte Carlo simulations indicate that = 
4.9 ± 0.7 events of these backgrounds survive the above selection criteria, which will be 
subtracted from the observed number of events in the calculation of the branching fraction. 
Background remaining in the continuum data sample is negligible. 

Fig. 2 shows the confidence level distribution from the kinematic fitting. The consistency 
between data and MC is satisfactory except for the first few bins in both Figs. 2(a) and 
(b). Figure 2(a) is similar to the situation for J/il) (see Fig. 4). Taking into 

account the similarity between the figures and that we expect little background in Fig. 4, 
we conclude that the discrepancy is due to the simulation of the error matrix in the track 
fitting rather than some unknown background, and this is included in the systematic error 
analysis. 
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TABLE I: Numbers used to calculate B['ip{2S) 3(7r+7r )]. 


"stse 


N>>9 


^^{2S) 

±3.686/±3.650 

670 

71 

4.9 

1.4 X 10^ 

6.8% 

3.07 


B. Branching Fraction and Systematic Error 


To obtain the branching fraction of 3(7r’''7r“), we have to subtract from iVgggg 

the number of continuum events at 3.686GeV. This number is estimated as fVg g^gg x /, where 
/ is the normalization factor; 

r ^cont 

r _ -^3.686 X <^3.686 
r \y ^cont ’ 

>^3.650 X <^3.650 

where, is the Born order cross section of continuum process, which is s dependent and 
can be expressed in terms of the 3(7r’''7r“) form factor 

a-%s) = ^x\T{s)\\ ( 1 ) 


where a is the QED hne structure constant. Assuming JF(s) oc 1/s, we get: 


j ^3.650,^g 


-^s.eso 


X — 
^3.686^ 

X 0.94. 


The branching fraction of t/’(25') 


-3.650 

3 (7r’''7r“) can be calculated as 


B[iIj{ 2S) 3{7i~^7i ) = 


N: 


3.686 


N: 


3.650 


X f 


£i>{2S) X A'^'( 2 S) 

= (4.83 ±0.38 ±0.69) x 10-^ 


where the hrst error is statistical and the second is systematic. The values of the variables 
in the equation are listed in Table I. 

The systematic errors come mainly from the MDC tracking, the generator, the continuum 
subtraction, and the kinematic fit, as well as the statistics of the MC samples. The systematic 
error contributions are listed in Table II and explained below. The total systematic error is 
15%. 


1. The MDC tracking efficiency was measured using J/'ip AA and 

n^n~ J/'il), J/'ip ^ events. It is found that the MC simulation agrees with data 



TABLE II: Summary of systematic errors. 


Source 

systematic error(%) 

MC statistics 

0.7 

MDC tracking 

12 

Kinematic Fit 

4 

generator 

4 

continuum subtraction 

4 

total number of 'ip{2S) 

4 

Total 

15 


within 1-2% for each charged track. The largest difference is taken as a conservative 
estimation, and 12% is quoted as the systematic error on the tracking efficiency for 
the channel of interest. 

2. Fig. 3 shows the 7r’''7r“ mass distributions after applying all the selection criteria; a 

clear p signal is observed in both the t/’(2S') and the continuum samples. To estimate 
the efficiency difference between the data and MC simulation (pure phase space), we 
generate samples with different intermediate states in 'ip{2S) 3(7r’''7r“), e.g., p, u, 

etc., and hnd the differences from pure phase space are at the 4% and 5% level for 
'^(25') decays and continuum data, respectively, which are taken as the systematic 
errors from the generator. 

3. Above, we assumed oc 1/s when subtracting the continuum contribution. As¬ 
suming a different dependency, such as JF oc 1/s^, yields a difference of 4% and is 
regarded as the systematic error of the continuum subtraction. 

4. The systematic uncertainty from the kinematic £t is estimated to be around 4% by 
using a different MDC wire resolution simulation model [9]. This is in agreement 
with various studies using pure data samples which can be selected without using a 
kinematic fit [12]. 

Compared to the previous Mark-I result [2] of 3(7r+7r“)] = (1.5 ± 1.0) x 10“^, 

our measurement has much better precision and a considerably higher central value. 


9 




FIG. 3; The 7r'’“7r“ invariant mass distributions after applying all selection criteria. The upper plot 
is for '4^{2S) 3(7r'’'7r“), while the lower one is for e+e“ ^ 3(7r'’'7r“). The dots with error bars are 

data, while the histograms are MC simulated phase space distributions. All vr^vr” combinations 
are included. 

C. Form Factor of e+e” —^ 3(7r+7r~) 


Using our values of fVs.gso and C and Eq. 1, the form factor at ^/s = 3.650 GeV can be 
determined. Following the procedure in Ref. [1], assuming only one photon annihilation in 
r/’(2S') —3(7r’''7r“) decays, one can also determine the form factor at ^/s = 3.686 GeV, using 
[ 1 ] 


Borni^s) — 


dTTQ;^ 

3s 

3s 


X |J^(s)px [l + 27^R(s) + |R(s)|2] 

x|.F(s)px[l + |R(s)p], 


( 2 ) 


with 


B{s) 


3y^rgg/G! 

^-^P(2S) +*^b(2S)rt’ 


where M^{ 2 S) is the mass of 'ip{‘2.S) and Fee and Ft are the partial width to e+e“ and total 
width of 'ip{2S). In Eq. (2), the interference term is neglected since it is at the 1.3% level [1] 
and small compared with the experimental uncertainties. 
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Using the numbers listed in Table I, we obtain 

I = 3.650^)1 = 0.21 ±0.02, 

I= 3.686^)1 = 0.20 ±0.01. 

The form factors measured at the two energy points are consistent with each other within 
la, i.e., there is no strong evidence for a large contribution from other than the electromag¬ 
netic interaction in 'ip(2S) decays. 

Using the form factor measured at the peak, the branching fraction of 

3 (7r’''7r“) is determined to be: 

B[ip{2S) 3(7r+7r")] = (4.8 ± 0.6) x 10■^ 

in good agreement with the result obtained by subtracting the continuum contribution 
determined from the continuum data, but with slightly improved precision since an extra 
assumption is introduced. In the above calculations, |JU(s)| oc ^ is assumed; assuming 
|jF(s)| oc ^ results in a difference in B[ip{2S) 3(7r’''7r“)] less than 4%, which is taken as a 
systematic error in this branching fraction. 

IV. MEASUREMENT OF J/tp 2(7r+7r-) 

As has been shown in Fig. 1, 'ip{2S) 3(7r’''7r“) hnal states can be used to measure the 

branching fraction of J/'^ —> 2(7r+7r“) via 'ip{2S) 'K^n~ J/ip, J/ip ^ 2(7r’''7r“). 

In order to decrease the systematic error, we determine the branching fraction of J/ip 
2{ti^ti~) from a comparison of the following two processes as has been done in Ref. [12]: 

• 1p(p2S) ^ J /ip, J /ip ^ 2{'K~^T^~) (I) 

• ip{2S) ^ 7r^7i~ J/ip, J/ip ^ (II) 

The branching fraction is determined from 

ATObsj 

B[J/ip 2(7r+7r")] = ^ x B{J/ip /iV), (3) 

Ajj /Eli 

where is the number of observed events, and e is the efficiency. The branching fraction 
for the leptonic decay J/ip ^ is obtained from the PDG [13]. 
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FIG. 4; The confidence level distribution for kinematic fitting of ^{2S) Tr~^ tt~ J/'ip, J/^|> 

2{Tr~^TT~) events. The dots with error bars are data, while the blank histogram is MC simulated 
signal. 

A. Event Selection 

For process I, we use the same selection criteria as for 'ip{2S) 3(7r+7r“) except the J/ip 

and KgKg vetos are no longer used. To remove t/’(2S') — n^7i~ J/tp, J/ip K*K + c.c. — 
+ c.c. background, events are vetoed if the recoil mass of one 7r’''7r“ pair falls 
into (3.07,3.12) GeV/c^ and, at the same time, the mass of another 7r’''7r“ pair falls into 
(0.47,0.53) GeV/cT 

Fig. 4 shows conhdence level distributions for the kinematic htting of ipi^S) —> 
7r~^7r~J/ip, J/ip The agreement between data and MG simulation is very similar 

to that of direct t/’(2S') —3(7r’''7r“) in Fig. 2(a). Because of the clean J/ip signal seen below, 
we may conclude that there is nearly no background in this process and that the discrepancy 
at small conhdence level is due to the simulation of the error matrix in track htting, which 
will be taken into consideration in the systematic error. 

For process II, the selection criteria are similar to those in Ref. [10]. The two lower 
momentum tracks are assumed to be 7r’''7r“, while the two higher momentum tracks 
The recoil mass of 7r'''7r“ candidates must fall within (3.0, 3.2) GeV/c^, while the invariant 
mass of candidates must be within 250 MeV/c^ of the J/ip mass. For candidates, 
each track must have > 2, where is the number of muon identihcation (MUID) 
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FIG. 5: 7r+7r“ recoil mass distributions of V’(2<S') ^ 3(7r+7r“) candidates. The histogram is MC 
simulation (upper plot), while the dots with error bars are data (lower plot). The smooth curves 
show the best fits to the distributions as described in the text. All 7r'’'7r“ combinations are included, 
and the background mainly comes from the incorrect combinations. 

layers with matched hits and ranges from 0 to 3, indicating not a muon(O), a weakly(l), 
moderately(2), or strongly(3) identified muon [11]. 

The recoil mass distribution of all tt+tt” combinations of t/’(2S') —3(7r’''7r“) candidate 
events is shown in Fig. 5. Both data and MC simulation are htted using a double-Gaussian 
for the J/'ip signal and a second order polynomial for the background. 

Fig. 6 shows the recoil mass distribution of the 'ipiplS) jj,'^iJ,~ candidates. 

A similar £t is performed as for the 3(7r’'“7r“) mode. 

From the fits we get = 0.027 ± 0.001 and ej/en = 0.46 ± 0.01. 

B. Branching Fraction and Systematic Error 

Using Eq. 3, the relative branching fraction of J/ip —>■ 2(7r’''7r“) is determined to be: 

^ = (5.86 ± 0.19 ± 0.39)%, 

F J 

where the first error is statistical and the second systematic. 

Using the branching fraction of J/'^ —>■ from the PDG [13], we obtain 

B[J/p^ 2(7r+7r-)] = (3.45 ± 0.12 ± 0.24) x 10"^ 
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Recoil Mass (GeV/c^) 

FIG. 6: The tt+tt” recoil mass distributions of '4>{2S) 7r'''7r“;U+^“ candidates. The histogram 

is MC simulation (upper plot), while the dots with error bars are data (lower plot). The smooth 
curves show the best fits to the distributions as described in the text. 

TABLE III: Summary of systematic errors for ^ 2(7r'’“7r“)]. 


Source 

systematic errors (%) 

MC Statistics 

1 

MDC tracking 

4 

Generator 

3 

Signal fitting 

1 

kinematic fitting 

4 

Muon identification 

2 

Background 

1 

Total error 

7 


Since this is a relative measurement, many systematic errors cancel, either completely or 
partially. The remaining systematic errors are listed in Table III and are described below. 

1. For process I, the statistical error of the MC sample is 0.9%, while for process II, it is 
0.7%, including any uncertainties introduced by the hts. 

2. There are two more tracks in J/ijj than in J/'ip fhe uncertainty 

in tracking is dominated by the two extra tracks and is estimated to be about 4%. 
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FIG. 7: The invariant mass distributions for J/'0 ^ 2(7r'’“7r“) candidates after applying all 

the selection criteria. The dots with error bars are data, while the histogram is MC simulated 
phase space. The dip at the left side of the p peak is due to the mass cut. There are four 
entries for each J/'0 ^ 2(7r'*'7r“) event. 

3. Fig. 7 shows the 7r’''7r“ mass distributions after all the cuts, where clear p and /2(1270) 
signals can be seen. To estimate the efficiency difference between data and MC sim¬ 
ulation (pure phase space), we generate samples with different intermediate states in 
J/ip 2(7r’''7r“), e.g., p, /2(1270), 02(1320), etc. and find the differences from the 
pure phase space MC simulation are about 3%, which is taken as the systematic error 
for the generator. 

4. The numbers of events obtained from the fits are affected by the background shape, as 
well as the signal shape. Using a linear background yields a change of the branching 
fraction of 1%, while using direct counting methods to estimate the numbers of events 
and efficiency results in negligible change in the branching fraction. The systematic 
error due to the fitting procedure is taken as 1%. 

5. The systematic uncertainty from the kinematic £t is estimated the same as that of 

'ip{2S) 3(7r+7r“) and is about 4%. 

6 . Instead of using the MUID, a sample may be obtained by using the energy 
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deposited in the shower counter to remove J/ijj e^e~ events and neglecting 
J/ijj —> TT+TT” and K^K~ which have small branching ratios [12], From the num¬ 
ber of events obtained with this method, we find that the error associated with the 
MUID requirement is about 2%. 

7. Background contamination is determined using Monte Carlo simulation. The back¬ 
ground fractions in both processes are far less than 1%, and the uncertainty of the 
estimation is at 0.5% level. 

Adding the systematic errors in quadrature, the total systematic error for B[Jl'ip —>■ 
2 (7r+7r“)] is 7%. 

Our measurement of B[J/ip 2(7r’''7r“)] is in good agreement with the existing measure¬ 
ment by Mark-I [3], which is (4.0 ± 1.0) x 10“^, based on 76 observed events. Although non 
resonance continuum 2(7r’''7r“) production was mentioned in the Mark-I paper, it is not clear 
whether it was removed in evaluating the J/ijj decay branching fraction. Our result has bet¬ 
ter precision compared with the Mark-I result [3]. A recent measurement of this branching 
fraction using initial state radiation events was reported by the BABAR experiment [4]; the 
result ((3.61 ±0.26 ±0.26) x 10“^) agrees with our measurement well but has a larger error. 


V. SUMMARY 


Using 14 M 'ip{2S) events and 6.42 pb“^ of continuum data at ^/s = 3.650 GeV, the 
branching fractions of 3{n^7i~) and J/ip are determined to be (4.83 ± 

0.38 ± 0.69) X 10“"^ and (3.45 ± 0.12 ± 0.24) x 10“^, respectively. The former is larger 
than the Mark-I result using 'ip{2S) decays [2], while the latter is consistent with Mark-I 
measurement [3] using a, J/ijj decay sample. In both cases, our measurements have improved 
precision. 

Using the above results and B['il){2S) 2(7r’''7r“)], B[J/ip 3(7r’''7r“)] from PDG [13], 

we obtain 


_ B[ip{2S) 3(7r^7r“) 

B[J/'iIj —> 3(7r^7r“)] 
_ B[iIj{2S) 2(7r’''7r“)] 

B[J/'ip 2(7r’'"7r“)] 


(12 ± 7)%, 
(13 ± 3)%. 


They are consistent with the “12% rule” [14] expectation within errors. 
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